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\q ', We obtain the most general static intersecting brane solutions by 

■^ | directly solving the relevant equations of motion analytically and 

in complete generality. These solutions reduce to the known ones 
in special cases, and contain further a class of new static solutions 
which are horizonless. We describe their properties and discuss 
their physical relevance. Along the way, we also describe the 
features of the brane energy momentum tensors, the equations 
of motion, and their solutions which arise as consequences of the 
intersection rules and the U duality symmetries of M theory. 
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1. Introduction 

Consider a static spherically symmetric star. If it is sufficiently massive 
then it is likely to collapse and to finally form a black hole or, possibly, some 
other static object. It is self evident that, generically, the fields outside such 
stars and the final collapsed objects must be described by the most general 
static solutions to the relevant equations of motion. In Einstein's general 
theory of relativity, for example, the most general static spherically sym- 
metric solution is given by the Schwarzschild solution. In Brans-Dicke the- 
ory, the most general static spherically symmetric solutions are given by the 
Janis-Newman-Winicour-Wyman (JNWW) solutions [1], the Schwarzschild 
solution now being a special case. 

In M theory (or equivalently string theory), the black holes are described 
by various intersecting brane configurations where the branes wrap the com- 
pact internal spaces, taken to be toroidal. For example, the four dimensional 
black holes are described in M theory by two stacks of Ml branes and two 
stacks of M5 branes which wrap the seven dimensional compact toroidal 
space and intersect according to the BPS rules. 1 The standard intersecting 
brane solutions are given in [2, 3, 4, 5, 6, 7]. They are static solutions in- 
dependent of compact toroidal coordinates and are spherically symmetric in 
the non compact transverse space. 

It is natural to expect that such black holes in M theory must have been 
formed by the collapse of sufficiently massive stars. In order to study such 
collapses, one first needs to construct a static star in M theory. This requires, 
among other things, a knowledge of the most general static solutions to the 
relevant equations of motion. 

The standard brane solutions given in [2, 3, 4, 5, 6, 7] are not the most gen- 
eral ones. These solutions can be thought of as analogous to the Schwarzschild 
solutions rather than to the JNWW solutions. Brane solutions of the JNWW 
type also exist. They have been found by several groups in different forms 
and using different methods [8, 9], see [10] also. However, in all these works, 
while solving the equations of motion, some ansatz or the other is made re- 



1 According to the BPS rules, two stacks of five branes intersect along three common 
spatial directions; a stack each of two branes and five branes intersect along one common 
spatial direction; and two stacks of two branes intersect along zero common spatial direc- 
tion. Also, the branes are taken to be uniformly smeared along the remaining compact 
directions [2, 3]. 



garding the form of the brane solutions which thereby limits their generality. 
Therefore, these solutions are not the most general ones. 

In this paper, we find the most general static intersecting brane solutions. 
We directly solve the relevant equations of motion analytically and in com- 
plete generality. No ansatz is made regarding the form of the solutions. The 
resulting solutions are thus completely general. In special cases, they reduce 
to the standard solutions of [2, 3, 4, 5, 6, 7] and to the more general ones 
of [8, 9]. The general solutions we find contain further a class of new static 
solutions which are horizonless and, to the best of our knowledge, have not 
appeared in the literature. 2 

We point out here that, except the Schwarzschild and the standard brane 
solutions, the general solutions mentioned above - namely, JNWW solutions, 
those given in [8, 9], and the new ones found in this paper - all have, gener- 
ically, a singularity either at a finite radius or at a vanishing radius. These 
singularities are not covered by any horizon and, hence, are naked. Also, 
they are present even though the solutions all have positive ADM mass. In 
the case of a static star, the interior solutions get modified and, hence, these 
singularities are not relevant since the general solutions are applicable only 
outside the star. The general solutions can then be used to explore the ob- 
servational consequences by studying the motions of various probes outside 
such a static star, as in [12, 13, 14, 15] for example. 

Also, such general solutions are essential ingredients in the study of col- 
lapse of stars, and in determining the nature and the properties of the final 
collapsed objects; in particular, in determining whether a naked singularity 
appears or not during and after the collapse. In the context of four dimen- 
sional JNWW solutions, such collapse processes and various issues involving 
naked singularities have been extensively studied in [14, 16]. Similar studies 
are also needed in the present context. 

In this paper, we describe our general static intersecting brane solutions. 
We start with the eleven dimensional low energy effective action for the inter- 
secting branes of M theory, and write down the equations of motion. These 
equations, their solutions, and the energy momentum tensors for the branes 
all exhibit several features which arise as consequences of the BPS intersec- 



2 In our earlier reports [11], we obtained the most general vacuum solutions first and 
then used it to generate the brane solutions by suitable boosting and U duality operations 
[6]. The brane solutions thus generated turn out to be the same as the ones obtained in 
this paper by directly solving the equations of motion. 



tion rules and the underlying U duality symmetries of M theory We describe 
these features along our way to solving the equations of motion. We first solve 
the equations and obtain the general analytical solutions in terms of a vari- 
able r which is suggested naturally by the equations themselves. We then 
obtain the general relation between r and the standard radial coordinate r, 
which completes the solutions. All along, we take particular care in ensuring 
that no generality is lost at any stage. We then describe the properties of 
the resulting solutions. 

The paper is organised as follows. In section 2, we present the action and 
the metric which can describe the static intersecting branes. In section 3, we 
write down the equations of motion in several convenient forms and discuss 
the structure of the solutions. In section 4, we describe the features which 
are consequences of the intersection rules and the U duality symmetries of M 
theory. In section 5, we solve the equations of motion in complete generality 
in terms of r . In section 6, we obtain the general relation between r and the 
radial coordinate r . In section 7, we describe the properties of the solutions 
pointing out the special cases which lead to the solutions known before, and 
the general cases which lead to the new solutions. In section 8, we make 
several remarks which illustrate various features of the present solutions. In 
section 9, we conclude with a brief summary and by mentioning a few issues 
for further studies. 

Three Appendices contain some useful expressions and formulae. In Ap- 
pendix A, we give the non vanishing components of the Riemann tensor for 
the metric used here. In Appendix B, following the analysis of [17], we give 
the general expressions for the ADM mass. In Appendix C, for the sake 
of completeness, we give the general vacuum solutions obtained earlier [11] 
directly in terms of r . 

A Note to the reader : The paper is somewhat long and contains many 
technical details which we have presented in order to solve the equations 
analytically and, at the same time, to make the generality of our approach 
obvious. Much of these details may be familiar to the reader. Hence, for the 
sake of the reader's convenience, we have collected the main results and the 
expressions for the general solutions together in one place in the beginning 
of section 7 before proceeding to discuss their properties in sections 7 and 8. 
So those readers who are mainly interested in the solutions may go straight 
to these sections at a first reading. 



2. General Set Up 

We study the static solutions describing M stacks of M theory branes, 
intersecting according to the BPS rules [2, 3] whereby two stacks of five 
branes intersect along three common spatial directions; a stack each of two 
branes and five branes intersect along one common spatial direction; and two 
stacks of two branes intersect along zero common spatial direction. We start 
with the low energy effective action for M theory branes and an appropriate 
static ansatz for the metric. 

The relevant part of the eleven dimensional low energy effective action 
describing M stacks of M theory branes mentioned above may be written in 
standard notation as 

where / = 1, 2, • • • ,J\f , 



Mi-Af„ /+2 



FmvM„ +2 is the (nj + 2)— form field strength for the I stack of nj— branes, 
and nj = 2 or 5 . We write the summation indices explicitly since not all the 
sums in this paper are over repeated indices. The above form of the brane 
action can be shown to follow as a consequence of the BPS intersection rules, 
see [7, 10]. Setting k 2 = 1 , the resulting equations of motion are given by 

7Zmn — n 9mn 7t = T MN = 2_^ Tmn(i) (2) 

1 i 

E d M (V=9 F MM2 - M ^+ 2 ) = (3) 



a; 
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where Tmn is the total energy momentum tensor and 
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is the energy momentum tensor for the I th stack of branes. 

We take the spatial directions of the brane worldvolumes to be toroidal 
and assume necessary isometries. Let the spacetime coordinates be given 
by x M = (t, x l , r, 6 a ) where i = 1, 2, • • • , n c and a — 1, 2, • • • , m . The total 
spacetime dimension D = n c + m + 2 = 11 . The coordinates x % describe 
the compact, n c dimensional, toroidal space; and the radial coordinate r and 
the spherical coordinates 9 a describe the non compact, (m + 1) dimensional, 
transverse space. In the following, we will assume that m>2 and study the 
static solutions. The line element ds which can describe the static intersecting 
branes may now be written as 

ds 2 = g MN dx M dx N = -e 2X °dt 2 + J2 e 2A ' (dx*) 2 + e 2X dr 2 + e 2 °dQ 2 m (6) 

i 

where dfl m is the standard line element on an m dimensional unit sphere. The 
non vanishing components of the field strengths are F 0il ... in r = d r A(i) where 
ii,---,i ni and Am denote, respectively, the spatial worldvolume directions 
and the (nj + 1)— form gauge field for the I th stack of branes. The fields 
(A , X 1 , A, a) and A(i) are functions of r only. 

It can now be seen that the energy momentum tensors Tmn(i) given in 
equation (5) are all diagonal. We denote these diagonal elements as 

^V)' ^V)' ^rW' Ta a(I)) = (~Pli Pili n /' Pal) ( 7 ) 

where p a i = pen fo r all a . The total energy momentum tensor is now given 
by T M N = diag (-p, p h II, p a ) where 

P = J2pi , Pi = J2Pu , n = J2 H/ , Pa = Pq = J2pm ■ ( 8 ) 
III I 

From the expression for Tmn(i) in equation (5), it follows that pu = p\\i 
or p±i according to whether the x % direction is parallel or transverse to the 
worldvolume of the I th stack of branes, and further that 



-Pi = P\\i - -P±i = Hz - -pai = 7 F oh-i ni r F Wl -^ r . (9) 



1 
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Note that, in the studies of stars or of cosmo logical evolution, one often 
assumes that the total energy momentum tensor is the sum of M individual 
energy momentum tensors each of which is diagonal and seperately conserved; 
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and then supplements the equations of motion with equations of state for 
each individual I which determine the diagonal components T M M <^ in terms 
of pi, or equivalently 11/ . In our case here, these properties of the energy 
momentum tensors follow from an underlying action for branes and from the 
BPS intersection rules, and equation (9) provides the equations of state for 
each individual I . 



3. Equations of motion 

We first write down the equations of motion for the fields 11/ and 
(A , A*, A, a) using equations (2), (4), (6), (7), and (8) only. In particular, 
we will not use equations (5) or (9). The resulting equations are, there- 
fore, valid quite generally wherever the metric is given by equation (6) and 
where the total energy momentum tensor is the sum of M energy momentum 
tensors each of which is diagonal and seperately conserved. 

Let 

a=(0,i,a) , A a = (A , A*, A a ) , p a i = (Poi,Pu,Pai) (10) 

where p Q j = —pj , and \ a = a and p aI = pen for all a . Define 

A = ^A a = A + ^A i + m ( j (11) 

a i 

Ti = £:r M M(/) = n / + Ejw • (12) 

M a 

It is straighforward to calculate the Riemann tensor corresponding to the 
metric given by equation (6). We have listed its non vanishing components 
in Appendix A. Using these components and the above definitions, it follows 
after some algebra that equations (2), (4), (6), (7), and (8) give 

(n 7 ) r = -n 7 Ar + E^ A r (is) 

a 

A'-E( A ") 2 = 2]Tll / e 2A + m(m-i)e 2A - 2 ' T (14) 

a I 

\ a rr + (A r - A r .) \ a r = ]T {- PaI + -^) e 2A + 5 aa (m - 1) e 2A - 2CT (15) 
where the subscripts r denote r— derivatives. 
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Equation (15) now suggests a change of variable from r to r given by 3 

e A dr = e A dr . (16) 

Then, for any function X(r(r)), we have 

X T = e A - x X r , X rr = e 2(A - A) (X rr + (A r - A r ) X r ) (17) 

where the subscripts r denote r— derivatives. The line element for the (m+1) 
dimensional transverse space given in equation (6) now becomes 

e 2A dr 2 + e 2 CT ^ = e 2a ^2 X ^2 + ^ (lg) 

where the field x is defined by 

X = A-a = A° + X)A i + (m-l)o- . (19) 

i 

Equations (13) - (15), expressed in terms of r, become 

(n 7 ) r = -HjAr + YsPalK (20) 

a 

A 2 -]T(A?) 2 = 2^n 7 e 2A + m(m-l)e 2x (21) 

a I 

Kr = E (-Pai + ^) e 2A + ^ (m - 1) e 2 * . (22) 

Equation for the field x defined in equation (19) can now be obtained. It is 
given by 

Xrr = Y. ( n * + Pm) e2A + (™ - !) 2 e ' X • ( 23 ) 

The above equations may be written more compactly upon defining G a/ 3 and 
G al3 by 

G a0 = 1 - cW , G a/3 = ^^ - ^ , (24) 



3 We note in passing that this variable r and the solutions obtained in terms of r bear a 
striking similarity to those which appear in the 'H-FGK formalism' [18]. Upon substitution 
of the solution for x( T ): see later in the following, the form of the (m + 1) dimensional line 
element given in equation (18) becomes an example of such a similarity; the Mp brane 
solutions given in equation (80) below is another example. 



and upon denning h a i by 

h a i = Yl G <*f> [-P0i + r>\ ) = n/ + PaI • ( 25 ) 

Then we get 

(n 7 ) T = -2UjA T + J2h a i\^ (26) 

a 

Y,G a p\ a T K = 2]Tn /e 2A + m(m-l)e 2 * (27) 

a/3 I 

\ a TT = Y, G ^ h Pi e2A + $ aa {m-l) e 2x . (28) 

Xrr = £/^/e 2A + (m-l) 2 e 2 * . (29) 

We require the solutions to correspond to asymptotically flat spacetime 
in the limit r — y oo . Therefore, the above equations of motion are to be 
solved with the asymptotic condition that, in the limit r — > oo, the fields 

(e*W >e *; A/)) "+ (l,l,l,r; ^) (30) 

where Qj are the charges. The definitions of (A, x) an d equation (16) then 
imply that, in the limit r — > oo , 

e A ^r m , e^r" 1 " 1 , (r - r) -* i (31) 

v 7 (m - 1) r"^ 1 v 7 

where Tq is a constant. Thus r ^- r from below in the limit r — >■ oo . 
Furthermore, since j- is positive, it follows that r decreases from tq as r 
decreases from oo . Hence t < tq for r < oo and, incorporating this feature, 
equation (16) may be written as 

dr= / dre x ' A . 



The constant To can be set to zero with no loss of generality. But we will not 
do so in the following in order to avoid the seemingly spurious negative signs 
in expressions involving (r — r) . 
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Note that, for a metric of the form given in equation (6) and which 
corresponds to asymptotically flat spacetime, one can obtain the ADM mass. 
In Appendix B, following the analysis of [17], we have given the corresponding 
general expressions for the ADM mass. 

Linear equations of state 

To solve equations (26) - (29), one requires the equations of state which 
give, for example, p a i = (poi , Pu , Pai) as functions of 11/ . Consider the case 
where these equations of state are linear and are given by 

h aI = U I +p aI = z I a U I (32) 

where z : a are constants and z* a = z^ for all a . Equation (26) can now be 
solved to give 

n 7 = n 70 e~ 2A+l1 — ► h aI e 2A = z a (n /0 /) (33) 

where IIjo are constants and I 1 are given by 

I 1 = I! z a X a ■ (34) 

a 

Using these expressions, equations for X a , I 1 , and \ become 

A° r = Y. z<Xl ( n ^ e 1 ') + S aa (m - 1) e 2x (35) 

/ 

l'„ = ^^(n J0 e F )+4m(m-l)e 2 ^ (36) 

j 

Xrr = E 4 (n/o /) +(m- l) 2 e 2 * (37) 

where 

z ai = j- G ^ ^ , gu = J2 G^zi z J p . (38) 

One may now express X a in terms of I 1 and x as follows. Solve equations 
(36) and (37) as simultaneous equations to obtain IT/o e l and e 2x in terms 
of l T TT and Xtt ■ Substituting these expressions into equation (35) gives an 
equation of the form 

K T = X« T , X a = Y^l I + a a X 
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where af and a' a are constants. It then follows that the most general solution 
for X a is given by X a = X a + L a r + a" a where L a and a" a are integration con- 
stants. Substituting this expression for X a into the definitions I 1 = J2 a z i^ a 
and % = A — a will lead to J\f+ 1 consistency constraints each on L a and a" a . 
A complete solution thus follows once equations (36) and (37) are solved for 
/ 7 (r) and x( T ) ■ Finally, equation (27) will lead to one further constraint 
among various integration constants appearing in the complete solutions for 

(I 1 , X, A a ) . 

As pointed out earlier, the various forms of equations of motion above 
have been obtained without using the expression for Tmn(i) given in equa- 
tion (5) or (9). Hence, these equations of motion are all valid quite generally 
wherever the metric is given by equation (6) and the total energy momentum 
tensor is the sum of M energy momentum tensors, each of which is diago- 
nal and seperately conserved. Thus, for example, these equations may also 
be used to study stars made up of M theory branes; similarly, with a few 
straightforward modifications involving t and r, and pj and Hi, they may 
also be used to study the cosmological evolution of a universe made up of 
intersecting M theory branes [19, 20, 21, 22, 23]. 

4. M theory branes 

Whether or not the equations for / 7 (r) and x( r ) can De solved explicitly 
depends on the values of the constants z^ and the consequent structure of 
Q IJ . It does not seem to be possible to obtain explicit solutions for general 
arbitrary values of z^ . Hence we now specialise to M theory branes and first 
incorporate the U duality symmetries of M theory which lead to a relation 
among z^ and to an elegant form for Q IJ . 

U duality relation among (Uj,p a j) 

We use the U duality symmetries here the same way as in our earlier 
works [20, 21, 22]. We note that, for a given A/", different intersecting brane 
configurations of M theory can be related to each other by suitable U duality 
operations - namely, by suitable dimensional reduction and uplifting to and 
from type IIA string theory, and T and S dualities in type IIA/B string 
theories. For a metric of the form given in equation (6), such an operation 
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leads to relations among \ a which in turn, through their equations of motion, 
imply relations among (U,p a ) . Although only time dependent cases in early 
universe were studied in [20, 21, 22], the U duality details of these works 
carry over and are applicable to the static cases also with only a few minor 
changes. The relevant details may be found in [20] and in Appendix A of 
[22]. Hence, we present only the main results here without repeating the 
details. 

It can be shown that the relations among (U,p a ) for intersecting branes 
of M theory, obtained by applying U duality operations as described above, 
are all satisfied if the individual (11/, p a j) obey the relation 

P||/ = IT/ + poi + p±i + m (pen - p ±I ) (39) 

where p»j and p±i are the pressures along the directions that are parallel and 
transverse to the worldvolume of the I th stack of branes. The above relation 
is a consequence of U duality symmetries and, therefore, must always be 
valid independent of the details of the equations of state. We further take 
Pm — P±i which is natural since the sphere directions are transverse to the 
branes. 4 

It can further be argued [20] that the equations of state may be written 
in terms of one single function only. For example, they may be written as 

h aI = Uj+p aI = z a T({*}j) , a = {0,||,±,ft} (40) 

where {*}/ denote brane quantities such as the number and the nett charge 
of the branes in I th stack, the constant coefficients z a and the functional 
dependence of J 7 on the brane quantities {*} are same for all I , and 

zn = z± , z\\ = z + z ± (41) 

as follows from pen = p±i and from equation (39). Note that equations (40) 
may still lead to equations of state of the type considered in equation (32) 
if T depends on II only and if the function .F(n) is such that, for example, 
J-'(n) = u^Ii or u^Ii for low or high magnitudes of II . Equation (32) 



4 We mention here that, in the studies of early universe or of stars, the constituent 
matter components also satisfy T r r > 1 . = T a a , 7 ,, equivalently IT/ = p a i . Using poi = —pi 

and p a i = pm = p±i, equation (39) then becomes p\\i = —pi + 2p±j [19, 20, 21, 22]. 
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then follows, with z^ = ZaU 1 where u 1 = u^ or u^ according to whether 
the magnitude of 11/ is low or high. 

Consider the case where z^ = z a u l and z a obey equations (41). It then 
follows from equations (38) and from BPS intersection rules that 



«i I {ni+l) _ . . \ . , 



= I 9 z + z ± - zA u l (42) 

g IJ = 2z [z± - z S IJ ) wV . (43) 

The form of z aI and the relations among them are consequences of U duality 
symmetries of M theory. The relations among z aI will be reflected in the 
solutions for X a when written in terms of I 1 , see equation (35) and subsequent 
comments. Note the elegant structure of Q IJ and its independence on m . 
These features of Q IJ are consequences of both the U duality symmetries 
and the BPS intersection rules. Similar Q IJ was also present in our time 
dependent early universe studies. 



Equations of state from action S 

Now consider the equations of state for M theory branes which are given 
in equation (9). They follow from the energy momentum tensor in equation 
(5), which in turn follows from the action S given in equation (1). Setting 
Poi = —pi, it follows from equation (9) that h a i = Hi + p a i are given by 

(hoi, h ¥ , h ±I , hm) = (2, 2, 0, 0) Ui . 

Clearly, these h a i satisfy the U duality relations in equation (39). They are 
of the form given in equation (40) with ^(j*}/) = 11/ and with z a given by 

(zq, z h z ± , z n ) = (2, 2, 0, 0) (44) 

which satisfy equations (41). Writing h a i = z^ 11/ with z^ = z^ 1 and 
u 1 = 1 , it follows from equations (42) and (43) that z aI and Q IJ for M 
theory branes are given by 

z 01 = zV = 2 {UI ~ 8) , z" = z™ = 2 {UI + 1} , Q" = -8 5" . (45) 
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Thus (*ll, z 1 -) = | (-2, 1) for Ml branes and (*», z x ) = | (-1, 2) for M5 
branes. 

5. Solutions in terms of r 

We now solve the equations of motion in complete generality for the case 
where the equations of state for M theory branes are given by equations (9). 
The fact that z^ = z± = in this case leads to crucial simplifications which 
in turn lead to explicit analytical solutions. 

Expressions for X a and A(i) in terms of I 1 and x 

Using (4,2{,4,4) = (2,2,0,0) and Q 1J = -8 5 IJ in equations (34), 
(36), and (37), we obtain 

I 1 = 2 (A + J2 X ') ( 46 ) 

ie\\i 

where % e ||j denote spatial directions of the I th stack of branes; and 

^ T = 8p 70 / , X rr = (m-l) 2 e^ (47) 

where pio = —Hio . These simplified equations for I 1 and % can now be 
solved explicitly. Integrating the above equations once, we get 

(lif - 16 PI0 / = 2 c 1 , { Xr f - (m - l) 2 e^ = | (48) 

where c 1 and c x are integration constants. Furthermore, it follows from the 
asymptotic conditions given in equations (30) and (31) that, in the limit 

e 1 ' -+ 1 , e* -> -^ . (49) 

(m - 1) (r - r) 

Using equations (47) for e l and e 2x , equation (35) for X a can now be 
integrated to give 

1 _ !) aa 

X a = --J2 zaI l I + -X + L a (r-r ) (50) 

8 j m — 1 
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where L a are constants, L a = L n for all a, and the asymptotic conditions 
given in equation (30) have been incorporated. Substituting the above X a 
into the definitions of I 1 and x leads to M + 1 consistency constraints on L a 
given by 

J2ziL« = , Y.L a = L a . (51) 

a ex 

Substituting \ a into equation (27), and using equations (48) and (51), leads 
to one further constraint among (c 1 ,c x ,L a ) given by 

S^^ilX-EG*^- (52) 

Note that the constraint J2 a L a = L n implies that 

afS a a 

= (L°) 2 + E(^) 2 + (^-l) (L n ) 2 >0 . 

i 

Hence c x > if £/ c 7 > . 

An expression for the field strength Fd^...^ r can also be obtained. Using 
Piq = —II jo, equation (33) for IIj, and equation (46) for I 1 , it follows from 
equation (9) that 

-2A+1 1 _ 1 „-2X-l ! ( iTi \2 



P/ = -n 7 = p 70 e-^ = - e—'' F t 
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Oii---i ni r) j 



and hence that pjo > . Furthermore, since Fo^...^ r = d T A{i) , it follows 
that 
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e 2A - 2A (^A/)) =(^A/)j =4p/oe 2 ' J 

where the first equality follows from equation (17). Using the asymptotic 
form for Am given in equation (30), we then get 

e A ~ A d r Am = d T Am = -2 ^pT e 1 ' (53) 

and, taking Qj to be positive with no loss of generality, 

2^p^=(m-l)Qj . (54) 
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Solutions for I 1 , \i an d A(i) 
We now solve equations (47) for I 1 and x ■ Let -X"(r) satisfy the equation 

X TT = be x =>- (X r ) 2 - 26 e x = 2c 

where 6 is a given constant and c an integration constant. X may be thought 
of as the coordinate of a particle moving in a 'potential' (— be x ) with 'energy' 
c . The second equation above implies that if b is negative then c must be 
positive. Thus, there are three cases depending on the signs of b and c . The 
corresponding general solutions may be written as 

e X = 2 2l \ 2 , b = a\ c = 2 7 2 



e A = 2 ! 2 , 6 = a" , c = - 2 7 ^ 



e X = 2 2 \ 2 , b = -a\ c = 2 7 2 
o^ cosh A y 

where y = 7(r* — r), a is obtained from b, and 7 and r* are two integration 
constants. Note that, in the above three cases, X is unbounded above and 
below, bounded from below, or bounded from above respectively. 

In the following, we take a and 7 to be positive with no loss of generality. 
The constant b is given, hence it determines the sign of b and the value of 
a . The constants 7 and r* may be determined, for example, as follows. 
Let (x,v) be the initial values of (X,X T ) at r = r in . Then v 2 — 2be x = 1c 
determines the sign of c and the value of 7 . The solution for e x can now 
be selected from above according to the signs of b and c . This solution 
evaluated at r = Tj„ then determines yi n = 7(7-* — Tj n ), and thus r* . 

Now consider equations (47) and (48). Let p/ > and c 1 be positive 
for all / . (We will later comment on the case where one or more of the c 1 
may be negative.) Then c x > as follows from equation (52). Hence let 
c 1 = 2(7 7 ) 2 and c x = 2j3 2 , and take j 1 and j3 to be positive with no loss of 
generality. Then, comparing with the equations and solutions for X given 
above, it follows that the general solutions for I 1 and \ may be written as 

e 11 = -. {lI) \ 2 , 2// = 7 / (r/ + r -r) (55) 

4 p I0 smh A yj 
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eX = 7 7^ ~T ' Vx = H^o- r) (56) 

(m — 1) smh y x 

where Tj are integration constants, and we have incorporated the asymptotic 
condition given in equation (49) for e x in the limit t ^- t . It follows from 
equation (52) that (3, 7 7 , and L a must satisfy the constraint 

-^/? 2 = ^E(7 7 ) 2 -E^^ • (57) 

in i. t j a/3 

We now express e l in a different form that facilitates comparisons with 
the brane solutions as presented commonly in the literature. Since e l — > 1 
in the limit r — > r , we have 



2 sink (7 7 T/) = 1 

T 

which may be satisfied identically upon setting 



T~ e ~ v i > — T~ e ~ °i 

where Ci = cosh 0/ and Si = sink Gj . It then follows that 



Pio n c y i TI Ci 

— = C I S I , e 7 - — 



Using the above expressions in equation (55), we now obtain 



e 



A? 



, # 7 = C 7 2 - e 27/ ( r - r °) Sj . (58) 



Using the expression for e l given above, equation (53) can now be integrated 
to obtain Am ■ Noting that Hj — > 1 and A(i) — > in the limit r — > Tq , and 
using ^P- = Cj Sj , it follows from equations (53) and (58) that 

Am = C -£ (l - **<"» >) ; (59) 

and from equation (54) for the charges Qi that 

(m - 1) Q/ = 2 7 7 C, 5 7 . (60) 
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Also, using I 1 = 2 j 1 (r — To) — 2 InHi , the solution for X a given in equation 
(50) can be written in terms of Hj as 

al Xaa ( al I\ 

** = Y, — lnH I + ^— iX+ (l«-£^J(t-t ) . (61) 

We note here that the functions Hi and A(i) and the form of X a given 
above turn out to be the generalisations of the corresponding harmonic func- 
tions, the (rii + 1)— form gauge fields, and the metric components which 
appear in the standard intersecting brane solutions as given, for example, in 

[2, 3, 5, 6, 7]. 

6. Solutions in terms of r : obtaining r(r) 

We now have the complete general solutions for static intersecting M the- 
ory branes, obtained by directly solving the equations of motion in terms of 
the variable r . To express these solutions in terms of the more familiar ra- 
dial coordinate r, we need to obtain the function r(r) in complete generality 
using equation (16). We place an emphasis here on the generality of r(r) 
because of the following. General solutions in terms of r have been obtained 
by several groups for systems which, at different levels, are more general than 
ours [10]. But, while translating such solutions to the r variable, invariably 
some ansatz is made which thereby limits the generality of the final solutions 
given in terms of r . 

Consider the line element ds given in equation (6). With no loss of gen- 
erality, we exchange the function e A for another function / by taking 

e 2X dr 2 + e 2 °dn 2 m 

thus e°"~ A = ta/7 • Using e A dr = e A dr and e A ~ a = e x , it now follows that 

dr 




r T = ^- = e A ~ A = rJf e x . (62) 

dr v 

Since e x — ¥ r m ~ l in the limit r — > tq , see equation (31), the general relation 
between r and r can be expressed as 

e x(r) = r m-\ e h(r) ( 63 ) 



where the function h{r) — > in the limit r — > tq . Using equation (56) for 
e x , it now follows that r{r) is given by 

r- 1 = e*- h = -J— —— fl . (64) 

to — 1 smh p {tq — t) 

Calculating r T from this expression and substituting it into equation (62) 

gives 

i — t h 

' f = — e~ x = cosh (3 (r — r) — sin/i /3 (r — r) . (65) 

' P 

Thus r(r) and /(r) are determined once /i(r) is known. 

Obtaining h(r) 

We now obtain h{r) first for vacuum solutions, thereby also illustrating 
our method in a simpler context; and then for brane solutions. The function 
h(r) turns out to be the same in both of these cases. 

Consider first the vacuum equations of motion in terms of r, namely 
equations (14) and (15) with 11/ = p aI = . Equation (15) gives 

A*' = b 1 ' F 

where i' = (0,i), A* = (A , A*), b % = (6°,6 J ) are constants, and the function 
F is the solution to the homogeneous equation 

F„ + (A r - A r ) F r = . (66) 

This equation for F implies that 

e A - x F r = F T = M =► F = M (t - r ) (67) 

where Ai = (to — 1) r™" 1 is an integration constant and we have incorporated 
the condition that F — > in the limit r — > r , see equations (30) and (31). 
As for a, note that there is a freedom in defining the r coordinate. Using 
this freedom, we set 

a = b n F + lnr 

with no loss of generality. 5 Thus, for vacuum solutions, we have 

\ a = b a F+ 5 aa In r (68) 



5 Upon substitution of cr(r), equations (14) and (15) will now give an equation for the 
function f(r) . See Appendix C. 
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where b a = (b°,b l ,b a ), and X a = a and b a = b n for all a . One can further 
set b Q = with no loss of generality, so that r becomes the conventional 
radial coordinate. However, we keep b Q non vanishing since the underlying 
structure is then clearer. Also, define 

B = Y,b a , K = -Y J G aP b a b p . (69) 

a a/3 

Using X a given in equations (68), we now get 

x = A- a =(m-l)lnr + (B -b n ) F . 

Comparing with equation (63) then gives h = (B—b n )F for vacuum solutions. 
Consider now the brane equations of motion in terms of r, namely equa- 
tions (13) - (15) with 11/ and p a i given by equations (9). Following the 
comments made below equation (38) regarding expressing X a in terms of I 1 , 
and following the analysis of section 5, it can be seen that A* = (A , A*) are 
now given in terms of I 1 by equation (50), or equivalently in terms of Hi by 
equation (61), where z aI are given in equation (45) and (r — To) is now to be 
replaced by the function F satisfying the homogeneous equation (66). Thus, 
we write 

X' = Y In Hi + b 1 ' F . 

i 4 
As for a, note that p a i = pen = P±i for the branes. It therefore follows that, 
using the freedom in defining the r coordinate, we can set 

z ±T 
a = J2 —T ln Hl + b n F + lnr 

with no loss of generality. Thus, for brane solutions, we have 

X a = Y^—lnHi + b a F+5 aa In r . (70) 

Using A a given above, we again get 

x = A-a=(m-l)lnr + (B-b n )F 
because it follows from equation (45) that 

J2 zaI - z m = (wj + 1) z¥ + (8 - rii) z A 
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~±i 



and, hence, the coefficient of the InHj term in x vanishes. Comparing with 
equation (63) then gives h = (B — b )F for brane solutions also. 
Writing F in terms of r using equation (67), we now have 

h(r) = k (r - T ) , k = M (B - b n ) (71) 

for both vacuum and brane solutions. The functions r(r) and f(r) then 
follow from equations (64) and (65), and are given by 

a p k (r -r) 

r m - 1 = -Jl (72) 

m — 1 sinh (3 (r — r) 

and 

/— k 

\J f = cosh /3 (ro — t) — — sinh f3 {tq — r) . (73) 

P 

Note that, in principle, one can now obtain the functions F and /, and 
thereby all other functions also, in terms of radial coordinate r by using 
equations (67), (72), and (73). However, it turns out that the functions F(r) 
and f(r) can be obtained more directly. It is straightforward to obtain their 
equations of motion in terms of r . We have earlier solved these equations 
analytically and in complete generality, and reported the results in [11]. For 
the sake of completeness here, we briefly describe these analytical solutions 
and their derivation in Appendix C. 

Relations between (b a ,A4) and (Z/*,7 7 ) 

The brane solutions given in terms of r are parametrised by the set of 
constants (L a ,'~f I ) , whereas those given in terms of r are parametrised by 
the set (b a , Ai) . We now obtain the relations between these two sets. 

Consider the expression for X a given in equation (70). Using equation 
(63) to write In r = ^j in this expression, substituting F = Ai (t — t ) and 
h — k (r — To) into it, and then comparing it with equation (61), we get 

sua ai „I 

Mb a - — k = L a -T\ —^ ■ (74) 

m- 1 j 4 y ' 

Using the expressions for z^ and z aI given in equations (44) and (45), the 
J\f + 1 constraints on L a given in equations (51) can be shown to imply that 
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k = M. (B — b ), which is already known, and that 

/ = M w 1 , w 1 = ]- J2 z* a b a = b° + ]T V (75) 

where i e ||j denotes the spatial directions of the I th stack of branes. Fur- 
thermore, using the b a given in equation (74), it can be shown that K defined 
in equation (69) is now given by 

-M 2 k = \ E(V) 2 - E^ l° rf - -^- k 2 • 

L I a/3 " i l 

Then the constraint given in equation (57) becomes 

£2 = k 2 + ^Zl m 2 K . (76) 

m 

Note that equation (75) gives 7 7 (r — r ) = w 1 F . Equations (58) and 
(59) then give e l = e ^ and 



*i 



72 „2w J F c2 /I . ^ ^ /'i „2w J F 



H I = C]-e^^S] , A/) = r ^ i (1-e^'J . (77) 

Using 7 7 = (m — 1) r™ _1 «; 7 , equation (60) for the charges Qi gives 

Q, = (2 w 1 Ci Sj) r™" 1 . (78) 

Note also that there is a scaling freedom, namely 

(F, M, b a ) — > (cF, c.M, c^ 1 fe Q ) 

where c is a constant, under which (6 Q F, .M6 a ) and hence the solutions 
(A Q , ii/j, ^4(i)) remain unchanged. Using this scaling freedom, we can set, 
for example, 2 (B — b n ) = 1 with no loss of generality. 

We now count the nett number of constants among each of the sets 
(L a ,7 7 ) and (b a ,A4) . The constants (L a ,7 7 ) are (n c + 2 + A/") in num- 
ber, and must obey (A/" + 1) constraints given in equations (51). Hence, 
there is a nett total of (n c + 1) constants among (L a , 7 7 ) . Similarly, (b a , Ai) 
are (n c + 2 + 1) in number but with no constraints. However, one can set 
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b Q = using the freedom in defining the r coordinate; and set, for example, 
2 (B — b n ) = 1 using the scaling freedom. Hence, there is a nett total of 
(n c + 1) constants among (b a , Ai) also. This equality of the nett totals is a 
further check that there has been no loss of generality in our derivation of 
r(r) . 

7. Properties of the solutions 

We have now completed the derivation of the general solutions for static 
intersecting M theory branes, expressed in terms of r or r . Below, for the 
sake of the reader's convenience, we first collect the main results and the 
expressions for the general solutions together in one place; and then discuss 
the properties of the solutions. 

Main results in terms of r 

In terms of the variable r, the general intersecting brane solutions are 
given by 

i 

and Am where 

z al Sjaa 

\ a = Y,-T lnH i + -X + L a (r-r ) , 

j 4 m — 1 



(m - 1) sink (r - r) ' ' Y 4 



Hl = C] - e 2 ^ J (— o) S j , Am = %^ (l- <**<?-*>)) , 

and the constants z aI are given in equation (45). In the above expressions, 
(/3, 7 7 ) are taken to be positive, j3 is given by equation (57), and L a obey the 
constraints given in equation (51). The general Mp brane solutions follow 
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upon setting M = 1 . Thus, for example, the line element for the general 
Mp brane solutions is given by 



ds\ Mp) = H Al] ds 2 p+1 + H A± e 2a (e 2 * dr 2 + dn 2 9 _ p ) (80) 

where the / subscripts are omitted, 

ds 2 p+l = -e 2Z °(— o) dt 2 + j- e ^(T-r ) ( ^ )2 

e 2a = e i+ e 2L« (r-ro) j ^ = L a _ £^7 

and (p, AH, A- 1 ) = (2, ^, |) for Ml branes and (5, =±, §) for M5 branes. 
The general vacuum solutions follow upon taking Si = , thus Hj = 1 and 
•4m = , and formally setting / y I = in the expression for L a . 

Main results in terms of r 

In terms of the radial coordinate r, the general intersecting brane solutions 
are given by 

ds * = _ e 2A« rft 2 + j- e 2X \dx 1 ) 2 + e 2 ° (^- + dQ 2 m ) (81) 



and A(i) where 



and 



x a = J2 V ln Hi + ba F + saa lnr ' 



Hl = Cj - e 2w ' F Sj , Am = %^ (l - e 2w ' F 



The functions F and / in these expressions are not obtained in terms r 
Instead r, F, and / are all obtained in terms of r and are given by 



o k (t -t) 
„m-l P e 



r =7 TY ■ ^7777 7' ^ = ^(r-r ), (82) 

[m — 1) smn p (tq — r) 
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k 

f = cosh (3 (t — r) — — sinh (3 (r — r) , (83) 

using which one can obtain F(r) and f(r), see the comments below equation 
(73). In the above expressions, (j3, w 1 , fc) are taken to be positive, 

777 — 1 

M = (m-1) r™- 1 , /3 2 = A; 2 + - — - M 2 K , 

m 

and (w T ,k,K) are defined in equations (75), (71), and (69). The general 
Mp brane solutions follow upon setting M = 1 . Thus, for example, the line 
element for the general Mp brane solutions is given by 



-2 _ ttAW ,2 , tjA 1 ^ J2b a F I a < , 2,o2 



rfs J Mp) = jf *" ds 2 p+1 + ff^ e 26 ^ ( ^ f r 2 ^ 2 _ p ) (84) 

where the I subscripts are omitted, 

ds 2 p+1 = -e 2b ° F dt 2 + Y.e 2blF {da 



t=i 



and (p, M, A 1 ) = (2, =p, |) for M2 branes and (5, =±, §) for M5 branes. 
The general vacuum solutions follow upon taking Sj = , thus H j — 1 and 

A/) = ° • 

Behaviour of (F, /, r) and (if/, *4(-0) 

We now describe the behaviour of the functions (F, /, r) and (i//, .A(-O) 
as r decreases from To to — oo . The behaviour of X a then follows. It is clear 
from their expressions that these functions remain finite and vary smoothly 
in the open range — oo < r < r . Therefore, we will focus on their behaviours 
in the limits where r — > 7~o and where r — > — oo . In the following, we set 
b Q = and 2 (B — b Q ) = 1 with no loss of generality. Then 

5 = 6° + ^6 i = - , 2k = M , K=(b°) 2 + J2(b i ) 2 -- ■ 



Asymptotic limit r — >■ r 
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By construction, the limit r — > tq corresponds to asymptotically flat 
spacetime. In this limit, it follows from equation (31) that 



r m x -)■ — ->■ oo . 



(m- 1)(t -t) 

rn — 1 

Hence M. (tq — r) — > ^-rr . It then follows that 

/ -+ 1 - k (t q - r) => f -+ 1 







j-m— 1 



„.m— 1 



e^ -> 1 + .M (r - r ) ->■ 1 



#/ -> 1 + (2w'S?) ^T , ^ 



y,m— 1 



Qi 



y.m—\ " \ ) j.m—1 



where Qj = (2w i CiSi) r™ _1 , see equation (78). Thus, the general brane 
solutions here all behave the same way as the standard ones in the asymptotic 
limit where r — > oo . 

Functions e F , Hj, and Am 

We now describe the behaviour of the functions e F , Hj, and Am as r 
decreases from tq and approaches the limit r —¥ — oo . 

The behaviour of the function e F is straightforward to see. Since F = 
M.(t — tq), it follows that e F varies smoothly and monotonically as r de- 
creases, and that e F — > in the limit r — > — oo . The function e 2w F behaves 
the same way as e F since w 1 is positive. Therefore, it follows that the 
functions Hj and Am vary smoothly and monotonically as r decreases and, 
further, that Hj — > C] and Am, — > ^ in the limit r — >■ — oo . 

Functions r(r) and /(r) 

We now describe the behaviour of the functions t(t) and /(r) as r de- 
creases away from To and approaches the limit r — > — oo . Their behaviour 
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depends on whether K = 0orK>0orK<0. We describe these cases 
seperately. 

K = 

Note that j3 = k when K = . Then, equations (82) and (83) give 

r m-l = tZ .ff = p - k (to-t) 

(m - 1) (1 - e~ 2k ( T °" r )) ' yJ 

Hence, / — > and r — >• r as r — >■ — oo . Indeed, using 2fc = M. = 
(m—1) r™ _1 , it follows that 

2k (r-ro) _ / _ -, _ r cT 



/ = e 2 ^ T - ToJ = e 1 



y.m—1 



The standard intersecting brane solutions as given, for example, in [2, 3, 5, 6, 
7] all follow from this K = case upon further choosing 2b° = 1 and U = . 
One then has 2«/ = 1 for all /, see equation (75), and hence e 2w F = e F and 
Qi = (CjSi) r™ _1 . It then follows that H I and A(j) are given by 

H - 1 + ^^ g J - ^^ 



More general static intersecting brane solutions have been obtained earlier 
by several groups in different forms and using different methods [8]. They 
can all be shown [9] to follow for other choices of b a , but still with K = . 

m — 1 

The solutions in [9] were obtained by taking / = e — 1 ^-rr as an 

ansatz; K = was then found as a constraint. The present solutions are 
completely general. They are obtained without making any ansatz and, 
hence, are applicable when K ^ also. 

K>0 

Note that /3 > k when K > . Then, equation (82) implies that, in the 
limit r — > — oo, 

r m-l _^ 2 P e (k-/3) (to-t) _^ q 

m—1 

Consider equation (83). It implies that, as r decreases from tq to — oo, 
\/J begins to decrease from 1, but remains strictly positive since cosh(*) > 
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4 sinh{*) when j3 > k . It also implies that, in the limit r — > — oo, 




It therefore follows that /(r) decreases from 1; remains strictly positive; 
reaches a non zero, positive minimum; and then increases to oo as r decreases 
from ro to — oo , equivalently as r(r) decreases from oo to . It can be shown 
easily that / reaches its minimum f min = 1 — 4j at r/ m j n where r/ m j n is given 

by tank (5 (r - T fmin ) = | . 

Thus, for K > 0, we have that / remains strictly positive, Hi remain > 1 
and finite, and < e F < oo in the interval < r < oo . Therefore, it follows 
that the K > solutions have no horizons. These general, horizonless, K > 
solutions are new and, to the best of our knowledge, have not appeared in the 
literature. These are the ones referred to as 'a class of new static solutions' 
in the title of this paper. 

K<0 

Note that f3 < k when K < . Then, equation (82) implies that, in the 
limit r —¥ — oo, 

r m-l ^ JA e (*-fl (ro-r) ^ ^ 
TO — 1 

Since r m_1 — >• oo in the limit r — )■ To, it now follows that r{r) decreases from 
oo; reaches a non zero, positive minimum; and then increases to oo again as 
r decreases from tq to — oo . It can be seen from equation (62), or by a direct 
calculation, that r reaches its minimum r min at r rm i n where f{r rm i n ) = . 

Consider equation (83). It implies that, as r decreases from r towards 
— oo, v7 decreases and must reach zero at a finite value r rm i n given by 
tanh (3 (r — T rmin ) = t • The right hand side of equation (83) then becomes 
negative for r < r rmin . However, we do not fully understand the solution 
for r < r rmin . We mention this K < case here only for the sake of 
completeness. Since it is not clear to us how to continue the solution for 
r < T rm in, nor how to analyse the equations of motion for the corresponding 
parameter values, we will not discuss the K < case any further in this 
paper. 
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8. Remarks on the solutions 

We will now make several remarks, mostly in the context of vacuum 
solutions. Since the functions Hi and A(i) are finite for r < To, these remarks 
may easily be adapted for the general intersecting brane solutions also. 

Consider general static vacuum solutions. The corresponding line element 
is given by 

ds\ Mp) = -e 2b ° F dt 2 + £ e 2b ' F (dxr + e 2bnp l^- + r 2 dn 2 m J (85) 



where r, F, and / are given by equations (82) and (83). 

771.— 1 

(1) Note that / = e = 1 — ^-rr for K = , and that the standard black 
n c — brane solution is obtained upon further choosing 2b° = 1 and b l = b n = 
. It is well known [6] that one can start with this black brane solution, 
boost it along a compact direction to first generate a momentum charge, and 
then use the boosted black brane to generate various M = 1 brane solutions 
by suitable U duality operations - namely, by suitable dimensional reduction 
and uplifting, and T and S dualities. Further repeats of boosting and U 
duality operations will generate other M > 1 intersecting brane solutions. 
The standard static intersecting brane solutions of string theory can also be 
generated in this way. 

More general brane solutions have been obtained earlier by several groups 
in different forms and using different methods [8]. We have shown in [9] 
that these solutions can all be obtained by suitable boosting and U duality 
operations, where now the starting vacuum solution is as given in equation 
(85), and still with K = 0, but with otherwise arbitrary values of b a . 

It turns out that the general static intersecting brane solutions presented 
here can also be obtained by suitable boosting and U duality operations. 
The starting vacuum solution is as given in equation (85), but with arbitrary 
values of b a irrespective of whether K = or not. We have presented in [11] 
the solutions generated by this method, and it is easy to verify that they are 
same as the ones given in this paper which have been obtained by directly 
solving the equations of motion in complete generality. Also, as shown in 
[11], the expressions for (B — b ) and K remain invariant under U duality 
operations. Hence, their values remain unchanged for all the solutions related 
by U duality operations. 
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(2) The appearance of hyperbolic functions in the solutions for e l , and 
thus for Hi, can now be related to the underlying boost operations mentioned 
above. In section 5 here, the Sink functions arise because the integration 
constants c 1 in equation (48) are taken to be positive. If c 1 were negative 
then trigonometric functions would have appeared in the solutions. 

Such solutions involving trigonometric functions can also be generated 
by U duality methods. Instead of a boost along a compact x % direction, one 
now performs a rotation in the (x^x- 7 ) directions, and proceeds with other 
U duality operations which will generate 'tilted brane' solutions, see [24] for 
example. The underlying rotation will lead to the appearance of trigono- 
metric functions in the solutions. In the present context, they correspond to 
choosing one or more of the integration constants c 1 to be negative. 

(3) Upon compactifying on the n c dimensional torus described by the 
coordinates x l , one obtains an effective d = m + 2 dimensional non compact 
spacetime described by the coordinates x M = (t, r, 9 a ) . The corresponding 
effective action and the line element may be written, symbolically and in the 
standard notation, as 

s {d) ~ Jd d x^^e AC (n {d) + E(^) 2 + ■■) , a c = £a* 
~ / d d x y^ (n (d] + £(d,y) 2 + • • •) 



and 



dsf d) = g (d) ^ dx^ dx v = -e 2X dt 2 + e 2X dr 2 + e 2a dVL 2 m 



where g( d \ ^ and TZm) are the d dimensional metric and the corresponding 

2A C 

Ricci scalar in the 'physical frame', and g( d ) ^ = e d - 2 g( d ) ^ and 1Z( d ) are 
those in the 'Einstein frame'. The d dimensional action S^ now contains 
n c number of scalars descending from the sizes of the compact dimensions, 
and M number of 1— form gauge fields descending from the (nj + 1)— form 
gauge fields A(i) ■ The couplings of the scalars among themselves and with 
the gauge fields are not the most general possible ones, but are dictated by 
the eleven dimensional parent action S given in equation (1). 

Clearly, the solutions presented in this paper lead to the most general 
static spherically symmetric solutions of the above d dimensional system with 
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scalars and 1— form gauge fields. For example, Janis-Newman-Winicour- 
Wyman (JNWW) solutions [f] follow now as a special case. These JNWW 
solutions can be viewed as Einstein frame solutions for the four dimensional 
Brans-Dicke theory or, equivalently, as solutions for a five dimensional theory 
with one compact coordinate; and, can be shown to follow from the present 
vacuum solutions upon setting m — 2, n c — 1, b n ^ and K = 0, see [9] 
also. 

(4) For K > 0, we have that / remains strictly positive, Hj remain 
> 1 and finite, and < e F < oo in the interval < r < oo . Thus, the 
K > solutions have no horizons. We note here that similar horizonless 
solutions have also appeared recently in several works [25]. The set ups and 
the contexts of these works are totally different from the present ones, and 
the reasons for the similarities of the solutions are not clear to us. 

(5) The ADM mass Madm for the solutions can be obtained from the 
expressions given in Appendix B. It is easy to see that Madm is positive for 
all the general intersecting brane solutions if b a satisfy the inequality 



2 {mb n + $>*)< 



in 



m — 1 



Using 2(B — b Q ) = 1, the above inequality becomes b n < b° + 2 ^-1) • We 
will assume that b a are choosen, e.g. b n = 0, IB = 1, and b° > , such that 
ADM mass is positive. 

Consider now the curvature tensors. Note that the Riemann tensor com- 
ponents Rabcd, calculated in the local tangent frame with indices A, B,- ■ •, 
correspond to tidal forces. The non vanishing components of Rabcd for the 
metric in equation (6) are given in Appendix A. If any of these components 
diverges then we assume conservatively that there will generically be a cur- 
vature singularity. Now, by construction, the limit t ^ t$ corresponds to 
asymptotically flat spacetime. Therefore, the components of Rabcd will all 
vanish in this limit. Furthermore, the fields behave smoothly and remain 
finite for all finite values of r < To ; hence, so will the Rabcd components 
for these values of r . 

Consider the limit r — > — 00 . Let the radial coordinate r — > r$ in 
this limit, hence rs = r$ if K = and r$ = if K > . Using the 
expressions given in Appendix A, it can be seen after some algebra that, 
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for the general static intersecting brane solutions presented here, except for 
those with 26° = 1 and b % = b n = , atleast one of the components of Rabcd 
diverges in the limit r ^t- r$ ■ Therefore, for generic values of b a , there is 
likely to be a curvature singularity at rs . Note that this singularity is not 
covered by any horizon, is naked, and is present even though the ADM mass 
is positive. This feature, that naked singularity is present even though ADM 
mass is positive, is also exhibited, for example, by JNWW solutions, by those 
in [8, 9], and by the recent ones found in [25]. 

(6) We will list a few reasons as to why the general static solutions given 
in this paper are useful and are also likely to be physically relevant despite 
the singularities at rs mentioned above. The present solutions are obtained 
by solving the equations of motion in complete generality. Also, they have 
positive ADM mass. Therefore, it is physically reasonable to expect that 
these solutions describe the general static configuration of the fields outside 
a static star in a D = n c + m + 2 dimensional theory with compact n c 
dimensional toroidal space or, equivalently, in a d = m + 2 dimensional 
theory with n c number of scalars. 

Some of the reasons for the usefulness and the physical relevance of the 
present solutions are as follows. 

(i) In the studies of the static or dynamical properties of a star, one must 
assume a most general initial field configuration in the exterior of the star, 
which corresponds to the most general solutions to the equations of motion; 
otherwise, one is likely to miss some crucial properties. Any restriction of the 
parameter ranges must then come from physical criteria, such as stability of 
the resulting static configuration of a star or as an attractor basin during its 
dynamical evolution; or from requiring consistency with observations. 

(ii) Let r* be the equilibrium radius of a static star. If r* > rs, which 
is trivially valid for K > case since rs = 0, then the singularities at rs 
mentioned above are not relevant since the equations of motion are different 
inside the star, namely for r < r* , and the present solutions are not appli- 
cable there. For r < r*, one then has to solve a different set of equations 
of motion which includes the stellar material, and then match the interior 
solutions onto the most general exterior ones. Further stability studies of 
the resulting static configuration of the star may put an upper bound on its 
mass, but they may or may not further restrict the allowed values of b a . 

(iii) If a star is sufficiently massive then it is likely to collapse. Let us 
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assume that a static configuration is reached at the end of such a collapse, and 
that the r = rs surface is exposed. Then the present solutions are applicable 
for r > r$ and, generically, a naked singularity is likely to be present at rs ■ 

It is now possible that the collpase dynamics is such that, irrespective 
of the initial values of b a prior to the collapse, the final values after the 
collapse are always given by an 'attractor basin' where 26° = 1 and b l = 
b Q = so that the r = rs surface is a standard non singular horizon. This 
may be accomplished, for example, if all the non trivial field configurations 
are 'radiated away' during the collapse. A singularity, similar to the one 
mentioned above or a different dynamical one, may or may not occur and it 
may or not may become naked during the process of collapse itself. 

In the context of four dimensional JNWW solutions, such collapse pro- 
cesses and various issues involving naked singularities have been extensively 
studied in [14, 16]. Similar studies are needed in the present context also. 

(iv) Note that the present solutions are obtained from a low energy 
effective action involving only two derivative terms. Such an action is unlikely 
to be applicable in the strong curvature regime; then, one has to invoke the 
underlying fundamental theory, namely string/M theory here. It is then 
reasonable to assume that singularities will be resolved within such a theory. 
If this is the case then the present solutions are indeed physically relevant, 
but are applicable only until the strong curvature regime is reached. 

(v) On the other hand, and more practically, one may simply assume that 
our universe is (n c + 4) dimensional with compact n c dimensional toroidal 
space, or that it is four dimensional with n c number of scalars; and that 
the present vacuum solutions describe the general static configuration of 
the fields outside a static star. It is then possible to study the motions 
of various probes around such stars, their orbitals, et cetera and compare 
them with observations. Such studies may lead to novel phenomena or other 
observational consequences. See [15] for such studies in the context of higher 
dimensions and Dp branes, and [13, 14] in the context of four dimensional 
JNWW solutions. 

9. Conclusions 

We now give a brief summary and conclude by mentioning a few issues 
for further studies. In this paper, we obtained analytically the most general 
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static intersecting brane solutions. These solutions reduce to the known 
ones in special cases, and contain further a class of new static solutions. 
We described the properties of these general solutions and discussed their 
usefulness and physical relevance. Along the way, we also described the 
features of the brane energy momentum tensors, the equations of motion, 
and their solutions which arise as consequences of the BPS intersection rules 
and the U duality symmetries of M theory. 

We now mention a few issues which can be studied further. 

Assuming that the present general solutions describe the fields outside 
the stars in our universe, one may calculate the motions of probes around 
such stars and study the observational consequences, as in [12, 13, 14, 15] for 
example. 

As mentioned in the Introduction, it is natural to expect that the black 
holes described in M theory have been formed by the collapse of sufficiently 
massive stars. It is desireable to study such collapse processes. Towards this 
goal, one may first construct a static star in M theory where the present 
general solutions describe its exterior fields. 

Once the static solutions for stars in M theory are obtained, their stability 
properties can then be studied along the lines of [26]. It is likely that a 
sufficiently massive star will collapse. Such a collapse may also be studied. 
In these studies of collapse, the M theoretic cosmological solutions we had 
obtained earlier in [20, 21, 22] are likely to be useful 

Appendix A : Riemann Tensor components 

Consider the metric given in equation (6) where i = 1, 2, • • • , n c and the 
fields (A , X\ A, a) depend on r only. Let 

dQ 2 m = h ab d6 a d9 b 

where a = 1, 2, • • • , m and h a b depend on 6s only. Let e M A denote the D— bein 
fields and let 

D _ C M N P Q r> 

J^ABCD — e A e B e C e D n MNPQ 

where Rmnpq is the Riemann tensor. It follows after a straighforward algebra 
that the non vanishing components of Rabcd for the metric in equation (6) 
are given by 

Rri'rj' = —Si'ji e {^X 1 ^ + (A* — A r ) A* J 

34 



Rrarb — ~h a b 6 [O'rr + (<7 r — X r ) <7 r ) 

Ri'i'k'v = (Si'i'Sj'k' — Si'k'SjH') e ^X l r X J r J 

Ri'aj'b = —Si'j'hab 6 X r <J r 
Rabcd = 6 ° Pabcd(h) + {Jl a dhbc ~ h ac hbd) 6 O r 

where i! = (0,i), X 1 = (A , A 4 ), and p a bcd(h) is the Riemann tensor corre- 
sponding to the metric h ab . 

The Riemann tensor components in terms of the variable r defined in 
equation (16) may be obtained from the above expressions by replacing (r, A) 
with (r, A) . Thus, 



Rri'rf = Si'y e- 2A (A;' r + (A;'-A T )A;' y 

Rrarb = ~h a b ^ [Prr + (o"r ~ A r ) <7 T ) 
Ri'j'k'l' = (Si'l'fij'k' ~ bi'k'&j'l') e \K-Kj 

Ri'aj'b — —Si'j'hab & A^_(7 r 

Rabcd = e a p a bcd(h) + (h ad h bc — h ac h bd ) e a T . 

Appendix B : ADM Mass 

Consider the metric given in equation (6) where the fields (A°,A*,A, cr) 
depend on r only and obey the boundary conditions given in equation (30) 
which correspond to asymptotically flat spacetime. The ADM mass Madm 
for such a metric is given, following the analysis of [17], by 

\ i ' r— >oo 

where u m and V Uc are the volumes of the m dimensional unit sphere and the 
n c dimensional torus. Let the asymptotic behaviour of (A , A*, A, a) in the 
limit r — > oo be given by 

( e -°, e* «», ^) - 1 - (c», <f, <f , c») ^ 
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where (c°, c l , c r , c ) and ro are constants. Evaluating the above expression 
for the ADM mass then gives 

M ADM = ^-^ (r^ 1 Cm) (86) 

where 

C M = m(c n - c r ) - (m - 1) (mc n + ^> l ) . (87) 

i 

We now write down the constants (c°, c l , c r , c n ) and Cm for vacuum 
solutions, for p— brane solutions , and for intersecting brane solutions. 

Vacuum solutions 

For vacuum solutions, we have 

' 2ct\ / c 2b n F \ 

2A° „2\ l „2\ e \ _ I „26°F ^2b'F e „26 n F 1 



e , e , e , 



9 1 l ' ' -f 



From the behaviour of the solutions described in section 7 in the asymptotic 
limit r — > To, it follows that 

m-l „m-l 

/ ->■ 1 , e -» 1 - 

^ fpm— 1 ' j,m— 1 

in the limit r ->■ oo . Hence, (c°, c\ c n ) = 2 (6°, 6*, fe c ) and c r = (2b n - 1) . 
It then follows from equation (87) that 

C M (vac) = m - 2 (m - 1) (m6 n + ]T 6*) (88) 

i 

for vacuum solutions. In the standard case, we have 

2(6°, b\ b n ) = (l, 0, 0) , C M{vac) = m . 

p— branes 

For M (string) theory p— brane solutions, the metric (in Einstein frame) 
is described by 

/ 2oA / 2b n F \ 

e 2A °, e 2A \ e 2 \ € —\ = Ih a " e 2b ° F , H A "° rA± e 2b * F , H A± ?-—, H A± e 2 ^ 
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where H = 1 + (1 - e 2wF ) S 2 with w = b° + £ie|| 6* and 5 = SWi 0, and the 
exponents A" is for directions parallel to the brane and A 1 - is for transverse 
directions. See equation (81) with Af = 1 and with n c > p in general. In the 
limit r — ¥ oo , the functions / and e F are as given earlier, and hence 



m-l 

Thus (e 2A °, e 2A \ e 2A , *£) -»■ 1 - (c°, c 4 , c r , c n ) ^ in the limit r ^ oo 
where, for a p— brane along (x 1 , • • • , x p ) directions, 



c 



•° - 26°-2 W A"5 2 



c* = 26 i -2«;A ll 5 2 /or i = l,---,p 

= 2b i -2wA ± S 2 for i = p+l,---,n c 

c n = 2b n -2wA ± S 2 

c r = 2b a -l-2wA ± S 2 . 

Noting that c n — c r = 1, and defining Xh = pA^ + (m + n c — p)A ± , it then 
follows from equation (87) that 

C M (p) = m - 2 (m - 1) (mb n + J2 b *~ wX h ^ 2 ) • 

i 

For M theory branes, m+n c = 9 and (^1, A- 1 ) = (-§, |) or (-§, f ) for M2 
or M5 branes. (For string theory p— branes m + n c = 8 and it can be shown 
that (AH, A- 1 ) = (z=l, 2±I) in Einstein frame.) In all these cases, A^-A 1 = 

-1 and (m + n c ) A 1 =p+l . Hence, X H = p(A^ - A 1 ) + (m + n c )A^ = 1, 
and we have 

C M (p) = m - 2 (m - 1) (mb n + Y^V - w S 2 ) (89) 

i 

for p— brane solutions. In the standard case, we have 

?M (p) 

Intersecting branes 
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2(6°, 6\ 6 n ; «;) = (!, 0, 0; 1) , C M {p) = m + (m - 1)S 2 



For M sets of branes intersecting according to the BPS rules, the p— brane 
analysis given above applies with only a few straightforward modifications. 
Performing the analysis, it can be seen that the final expression may be 
obtained upon replacing wS 2 in the p— brane expressions by J2i w I S 2 where 
I = 1, 2, • • • , M and w 1 are given in equation (75). Thus, 

C M (jv P ) = m - 2 (m - 1) (mb n + £ V - £ «/ 5 7 2 ) (90) 

i J 

for intersecting brane solutions. In the standard case, we have 

2(6°, b\ b n ; w I ) = (l, 0, 0; 1) , C M {m = m + (m - 1) ^TS, 2 . 



Appendix C : Anayltical solutions for the functions F and / 

The equations of motion for F and / in terms of r may be obtained by 
substituting X a given in equation (68) into vacuum equations of motion (14) 
and (15) with 11/ = p aI = ; or, by substituting X a given in equation (70) 
into brane equations of motion (13) - (15). After some algebra, the resulting 
equations for F(r) and /(r) can be written as 



2(B - b n ) f (rF r ) -- 


- { m -l){l-f) + -f{rF r f 
m 


(91) 


2(B-b Q )f(rF r ) = 


= 2(m-l)(l-/)-r/ r 


(92) 


e (B- b -)F Jj {rFr) _ 


M 


(93) 



where Ai = (m— 1) r™ -1 and B and K are defined in equation (69). Equation 
(93) follows directly from equation (67), namely from e A_A F r = Ai . We 
now set 2(B — b ) = 1 with no loss of generality. For K = 0, the above 
equations can be solved easily and the solutions are given by 



m-l 

/ = e* = 1 - ' 



,-P _ -i r o 



j.m—1 



It turns out that equations (91) - (93) can be solved analytically for 
K > also; and, that the properties of the functions F(r) and f{r) can 
be understood even without obtaining solutions. These properties and the 
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explicit analytical forms of the functions F and / have been described in our 
earlier reports [11]. For the sake of completeness, we present these solutions 
here along with a description of their derivation and some of their properties. 
Define R = r m_1 and i? = r™ -1 . We then get 

(m- l)K 



f(RF R ) = 1-/+^ —f(RF R) 



m 
f(RF R ) = 2(1 -f)-Rf R 



e F f(RFn) 2 = § 
where the subscripts R denote R— derivatives. It follows that 

(1 _ ,) _ Rfs = (J^IE f (RFR f - (•» - «* W - /) - */«> 



.1 \, •—■ IX j. 

m m j 

which gives a quadratic equation for Rf R , and that 



e F 



Rlf 



it> 2 (2(1 - /) - Rf R Y 

Defining 

, 4(m-l)K 1 

b = , a = — — , /o = 1 - a 

m l + o 

and solving the quadratic equation for Rf R gives, after some algebra, 

/ " A + e \/a / (/ - /o) 

and 

,2 



6 4 a i? 2 (1 - /) 2 l j 

where e = ±1 and the square roots are always to be taken with a postive 
sign. 

We consider the case K > 0, hence 6 > 0, so that a < 1 and / > . 
Define R min by f{R min ) = /o , and a function #(i?) by 

y/f-fo = tgV® 9 => /=l-a + afi' 2 (96) 
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where e g = Sgn g and, further, choose e g = e with no loss of generality. 
Using the asymptotic behaviour of /, namely / — > 1 — -^ as i? — >■ oo , it 
can be seen from equation (94) that e = +1 , and from equation (96) that 
g(R) — > 1 — 2^ , i n this limit. It can further be seen that g > and e g = +1 
for R m i n < R < oo ; that g(R m in) = ; and that g < and e g = — 1 for 
< R < Rmin . Writing / in terms of g and after some algebra, equations 
(94) and (95) now become 

Rg R =l^L= =► ^ = d g[^^L) (97) 

and 

(98) 

It turns out that equation (97) can be solved and an explicit analytical solu- 
tion for R in terms of g can be obtained, because 

9 + ^/f = _J_ «_ 1 f 1 _M J_ 

l-<? 2 l-^ 2 v/J 2 ^l-<? 1 + J 77 

and each term on the right hand side can be integrated in a closed analytical 
form. Incorporating g(R m in) = , equivalently R(0) = R m in , we get 



1-2 2 

9 + v7 


dR , fg + y/J 

R =dg [i- 9 > 


f? 2 

r F _ n 

4a 2 R 2 


(g + VfV 



In 



l min 1 , M 2| i /~ o • L— 1 ^ V Q 



1 ^ 1 -« + «# 
— MWI ^=^ 

+± s^- 1 l ~ a ~ ag 



a) 



1 +g) Jail -a) 



A more explicit expression for R(g) can be obtained, and in many equivalent 
forms, by using identities such as 



Sink x x = In (x + vl + x 2 ) 



f-g 2 = (l-a)(l-g 



,2 



a (1 - a) (1 ± gf + (1 - a =p a#) 2 = / 
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(g-\ff) (i-«-^-v / 7) = (i-a)(i + <?)(i-v / 7) 

(l-a-ag + y/Tj (l-a + ag + y/Tj = (1- a) (l + \/T) 



The expression we find convenient is given, after a series of manipulations, 
by 

Rrain y/l - OL (1 - #) (1 + \fj) { Vf + 9V®\ / gg x 



-a (l-g)(l + y/f) f yff + gVa X 
1 - a + ag + y/J \ \/l — a J 



R 

Equation (98) now gives an expression for e F in terms of g . 

We now describe some properties of the analytical solutions given above. 
These properties all follow straightforwardly, but after some algebra. 

• The above expression for R(g) satisfies the differential equation (97). 



Consider a = 1 . Noting that ( y/T^~a) = 1 in the limit a — > 1 , 

and taking y/J = g , we get ^ = 1— g 2 = 1 — f . It then follows from 

v - f - 1 

R 



equations (97) and (98) that e F = f = 1 - ^ . 



Consider a < 1 . As R decreases from oo to R m i n to R± to , 5r de- 
creases from 1 to to —1 to — oo ; hence, / decreases from 1, reaches its 
minimum /o, then increases and reaches 1 again, and increases further 
to oo . Thus, R(g) for some select values of g are: 

fl(l_) = oo , R(0) = R mm , R(-l) = R 1 , R(-oo) = . 

As R decreases from oo to , e F decreases monotonically from 1 to 0, 
remaining < / always. 

Consider the limit R — > oo . Setting 1 — g = -^^ and f — g — 1 in 
equation (99) gives 

c (a) 1 / \ 1+ ,f* / \ ^Z" 

Rmin = ^0 , C(a) = - (1 + y/aj ~ [1 - y/aj 

which thus expresses Rmi n , where the function / reaches its minimum, 
in terms of Ro = r™ -1 . Note that \ < c(a) < 1 for < a < 1, as can 
be shown easily. 
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• Consider R = R\ . Then g = — 1 and / = 1 . Equation (99) now gives 

Ri = c(a) R min . 
which thus expresses R\ in terms of R m i n and thereby in terms of Rq . 

• Consider the limit R — > . Then g — > — oo , v7 ~~ * ( — 9)V® ) an d 
yj + ^^ = ^^a! ► i=« i . we then get 

i?- 1 ~ (-t/) 1 -^" , / ~ g 2 ~ i?"^ . 

• In the main body of the paper, the general solutions for R = r" 1 " 1 , F, 
and / have been obtained in terms of r , and are given in equations 
(82) and (83). Using equation (96), one can now obtain an expression 
for g also in terms of r . It can then be shown, after a long but 
straightforward algebra, that these expressions for R(t), F(t), f(r), 
and g(r) satisfy equations (98) and (99). 
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